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CN ■ Abstract 

(-h \ Thermodynamic properties of strange quark matter in strong magnetic fields H up to 10 20 G are considered 

within the MIT bag model at zero temperature implying the constraints of total baryon number conservation, 
charge neutrality and chemical equilibrium. The pressure anisotropy, exhibiting in the difference between 
the pressures along and perpendicular to the field direction, becomes essential at H > H t h, with the estimate 
10 17 < H t h 10 18 G. The longitudinal pressure vanishes in the critical field H c , which can be somewhat 
less or larger than 10 18 G, depending on the total baryon number density and bag pressure. As a result, 
£L[ the longitudinal instability occurs in strange quark matter, which precludes: (1) a significant drop in the 

content of s quarks, which, otherwise, could happen at H ~ 10 20 G; (2) the appearance of positrons in weak 
53"! processes in a narrow interval near H ~ 2 • 10 19 G (replacing electrons). The occurrence of the longitudinal 

■ instability leaves the possibility only for electrons to reach a fully polarized state, while for all quark flavors 
the polarization remains mild even for the fields near H c . The anisotropic equation of state is determined 

■ under the conditions relevant to the interiors of magnctars. 
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Strange quark matter, composed of deconfined u, d and s quarks, can be the true ground state 
of matter, as was suggested in Refs. [lrlsj]- There it was found that, at zero temperature and 
pressure, the energy per baryon in strange quark matter for a certain range of the model QCD- 
related parameters can be less than that for the most stable 56 Fe nucleus. This conjecture, if 
will be confirmed, would have important astrophysical implications. In particular, strange quark 
matter can form strange quark stars self-bound by strong interactions [J-lg]. This is in contrast 
to an ordinary scenario, in which neutron stars are composed of hadrons (plus some admixture of 
leptons to ensure charge neutrality and beta equilibrium) and are bound by gravitational forces. 
Also, if strange quark matter is metastable at zero pressure, it can appear in the high-density 
core of a neutron star as a result of the deconfinement phase transition. In this case, the stability 
of strange quark matter is provided by the gravitational pressure from the outer hadronic layers. 
Then a relevant astrophysical object is a hybrid star having a quark core and the crust of hadronic 
matter. Moreover, strange quark matter can be potentially encountered in the form of small 



nuggets called " strangelets" 



, Q] • For the review on the properties of strange quark matter and 
its possible forms in astrophysics, one can address to Refs. 

Another important aspect related to the physics of compact stars is that they are endowed with 
the magnetic field. In particular, for conventional neutron stars the magnetic field strength at the 
surface can reach the values in the range of 10 9 -10 13 G [lo| . For the special classes of neutron 
stars such as soft 7-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs), the field strength 
can reach even larger values of about 10 14 -10 15 G 111. Il2|. These strongly magnetized objects are 
called magnetars 13] . It was also suggested that a magnetized hybrid star, or a magnetized strange 



quark star can be a real source of the SGRs or AXPs 14|, LL5|]- The actual mechanism, by which 
magnetars generate such strong magnetic fields, is still under debate. Together with the dynamo 
amplification scenario in a magnetar with the fast rotating core (111 ], other possibilities such as 



spontaneous ordering of hadron Jig], or quark [l7|] spins are not excluded. 

Under such circumstances, the issue of interest is the impact of a strong magnetic field on ther- 
modynamic properties of neutron star matter 1814241]. or hybrid/quark star matter 2514291] . Note 
that in the interior of a magnetar magnetic field strength can reach even larger values compared 
to that at the surface. In the recent study it was shown that, for either of the scenarios 
with a gravitationally bound neutron star or with a self-bound quark star, the field strength in 
the magnetar core could be as large as 10 20 G. In such ultrastrong magnetic fields, the effects 
of the 0(3) rotational symmetry breaking by the magnetic field become important [liJQ. In 
particular, the longitudinal (along the magnetic field) pressure is less than the transverse (perpen- 
dicular to the magnetic field) pressure resulting in the appearance of the longitudinal instability 
of the star's matter if the magnetic field exceeds some critical value. The effects of the pressure 
anisotropy should be accounted for in the consistent study of structural and polarization properties 
of a strongly magnetized stellar object. As a consequence of these effects, the equation of state 
(EoS) of the stellar matter becomes essentially anisotropic in an ultrastrong magnetic field. 

In the given research, we consider the effects of the pressure anisotropy in strange quark matter 
under the presence of a strong magnetic field. To study quark matter, we use a phenomenological 
MIT bag model [34J], in which quarks of various flavors are considered as degenerate Fermi gases, 
confined in a finite region of space by the inward bag pressure. The MIT bag model provides 
simple and physically transparent framework to approach the problem and was applied earlier to 
the study of nonmagnetized 0, IE, 0] arid magnetized 25-29] strange quark matter. Note that the 
effects of the pressure anisotropy were not included in the research of magnetized strange quark 
matter in Refs. (2514281] . Besides, as we will discuss in more detail further, the determination of the 
total anisotropic pressure in Ref. |29J] was incomplete, where, in fact, the pure field contribution 
(the Maxwell term) was missed. However, just this term provides the most principal source of the 
pressure anisotropy in strong magnetic fields H > H t h, where 10 17 < H t h < 10 18 G (Hth is the 
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threshold field at which the pressure anisotropy becomes relevant). As a result, the longitudinal 
pi and transverse pt pressures interchange their roles: instead of inequality pt < p\ obtained in 
Ref. |29|], one gets the opposite, pt > pi- 

It is worthy to note at this point, that at sufficiently high density strange quark matter will 
be in color superconducting color-flavor-locked (CFL) state [35l. l36|. in which quarks of all flavors 
and colors near the Fermi surface are paired. Nevertheless, it is unknown which of the phases, 
color superconducting or normal, will be preferable at the densities closer to the density of the 
deconfinement phase transition. Recent perturbative QCD calculations up to the second order on 
the coupling constant a s with the finite strange quark mass still do not allow nor confirm nor rule 
out the existence of unpaired strange quark matter, if to take into account the uncertainties in the 
model parameters 37[. By this reason, the study of the impact of a strong magnetic field on the 
thermodynamic properties of normal strange quark matter is expedient. 



I. GENERAL FORMALISM 

In the simplest version of the MIT bag model, quarks are considered as free fermions moving 
inside a finite region of space called a " bag" . The effects of the confinement are accomplished by 
endowing the finite region with a constant energy per unit volume, the bag constant B. The bag 
constant B can be also interpreted as the inward pressure - the "bag pressure", needed to confine 
quarks inside the bag. In the MIT bag model, the bag constant is considered as a phenomenological 
parameter of the theory. Therefore, all relevant equations can be obtained, first, by considering 
the relativistic degenerate system of free fermions in an external magnetic field, and then, in order 
to get the equation of state of the system, by properly modifying the energy and the pressure of 
magnetized noninteracting fermions with account of the bag constant. 

The Lagrangian density for the relativistic system of noninteracting quarks (u, d and s) and 
leptons (e) in an external magnetic field reads 

£= Y, U^Wv-QiAJ-rrHtyi-^F^F^. (1) 

i=u,d,s,e 

Further we will assume that the external magnetic field is directed along the z axis, and, 
correspondingly, choose the 4-potential in the form = HxiS^z (/i = 0,1,2,3), where H is 
the magnetic field stren gth . In Eq. ([1]), = d^A u — dyA^ is the electromagnetic field tensor. 
Analogously to Refs. 2^-2^], we did not take into account in Eq. (pQ) the quark anomalous magnetic 



moments which are not completely understood in the deconfined state. Also, we disregarded the 
electron anomalous magnetic moment, whose role, as believed, is insignificant even for the magnetic 



fields encountered in magnetars [38]. The quark ipf (/ = u,d,s) and electron ip e spinors satisfy 
the Dirac equation 

[^{idfj, - qiA^) - mi]ipi = 0, i = u,d,s,e. (2) 
The energy spectrum of free relativistic fermions in an external magnetic field has the form 



I 1 s 

4 = \J k z + m i +2v\qi\H, i/ = n + ---sgn(^), (3) 

where v = 0,1,2... enumerates the Landau levels, n is the principal quantum number, s = +1 
corresponds to a fermion with spin up, and s = — 1 to a fermion with spin down. The lowest 
Landau level with v = is single degenerate and other levels with v > are double degenerate. 
For positively charged particles, the lowest Landau level is occupied by fermions with spin up, and 
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for negatively charged particles by fermions with spin down. As a result, each charged fermion 
subsystem acquires spin polarization in a magnetic field. 

Further we will consider thermodynamic properties of magnetized strange quark matter at zero 
temperature. This approximation is quite reasonable taking into account that for the characteristic 
densities in the interior of a neutron star of about several times nuclear saturation density the 
temperature is much less than the quark chemical potentials. Concerning electrons, although finite 
temperature effects can be important for them, their contribution to the thermodynamic quantities 
is usually unessential because the electron fraction, in turn, is small. In the zero temperature case, 
the thermodynamic potential for an ideal gas of relativistic fermions of ith species in the external 
magnetic field reads [iH] 



where the factor (2 — 5„o) takes into account the spin degeneracy of Landau levels, gi is the 
remaining degeneracy factor [gf = 3 for quarks (number of colors), and g e = 1 for electrons], /ij is 
the chemical potential, and 



111; 



mf + 2v\ qi \H , k F = Jul -m? . (5) 



In Eq. @, summation runs up to 



"max 



2\ gi \H J ' 



/[...] being an integer part of the value in the brackets. The number density Q\ = —{-qj^)t °f 
fermions of zth species is given by 



\li\mH 



Y(2-5 u , )k* (6) 



2vr2 ^ 

The sum in Eq. ([6]) can be split into two parts representing the fermion number densities with 
spin up and spin down. As explained earlier, the only difference between the two sums is in the 
term with v = 0, corresponding to spin-up fermions if they are positively charged, and to spin- 
down fermions, if they are negatively charged. Then the zero temperature expression for the spin 
polarization parameter of the ith species subsystem reads: 

Qi-Qi Qi9iH [~ 2 — 2 



n * = ^— "2^^-^- (7) 



In a strong enough magnetic field, when only a lowest Landau level is occupied by fermions of ith. 
species, a full polarization occurs with |Ilj| = 1. 

In order to find the chemical potentials of all fermion species (and, hence, the corresponding 
particle number densities), we will use the following conditions. First, the conservation of the total 
bar yon number is implied: 

^(Qu + Qd + Qs) = SB, (8) 



4 



where qb is the total baryon number density. Further, quark matter is considered to be charge 
neutral: 



2£»n - Qd~ Qs~ 30 e - = 0- 
Because of the weak processes in the quark core of a neutron star 0] 

d — > u + e~ + P e , u + e~ — > d + u e , 
s — > u + e~ + v e , u + e~ — > s + u e , 
s + u f-> d + u, 



(9) 



(10) 

(11) 
(12) 



all fermion species are assumed to be in a chemical equilibrium with the corresponding conditions 



M<2 = P-u + Me- ) 
= Ms- 



(13) 
(14) 



Here we suppose that neutrinos and antineutrinos freely escape a neutron star, and, hence, their 
chemical potentials are set to zero. Eqs. ©, ([9]), (|13p and (|14|) . with account of Eq. ([6]), form 
the full set of the self-consistency equations for finding the chemical potentials /ij of quarks and 
electrons. 

At zero temperature, the energy density E{ = fij + //jgj for fermions of ith species reads 



Ei 



4vr 2 



(15) 



177; 



In the MIT bag model, the total energy density E, longitudinal pi and transverse pt pressures 
in quark matter are given by 



E = Y,Ei + ^ + 



8vr 



Pi 



Pt 



^ 1 Stt 



H 1 



Vfii- iTM + B, 



(16) 
(17) 
(18) 



where B is the bag constant, and M = Yli M{ = ~~ Yli(7fit)iH i s the total magnetization. At zero 
temperature, the magnetization of ith fermion species reads 



Mi 



mm 

4vr 2 



(19) 



u=0 



x (A«i4> - + 2i/|gi|fl"] In ^ + ^ ). 



It is seen from Eqs. (I17j) . (118p that the magnetic field strength enters differently to the longi- 
tudinal and transverse pressures that reflects the breaking of the 0(3) rotational symmetry in a 
magnetic field. In a strong enough magnetic field, the quadratic on the magnetic field strength 
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term (the Maxwell term) will be dominating, leading to increasing the transverse pressure and to 
decreasing the longitudinal pressure. Hence, there exists a critical magnetic field H c , at which the 
longitudinal pressure vanishes, resulting in the longitudinal instability of strange quark matter. 
Further we will find the critical magnetic field H c for the appearance of the longitudinal instability 
in a quark system under the conditions relevant to the cores of magnetars. Also, we will deter- 
mine the threshold magnetic field beyond which the pressure anisotropy in strange quark matter 
becomes significant and cannot be disregarded anymore. 

Note here a principal difference between Eqs. (fTT[) . ([18]) . used in the present study, and analogous 
equations in Ref. [29(. Namely, the pure field contribution (the Maxwell term) was missed in 
the equations of Ref. and, by this reason, the authors of that work wrongly arrived at the 
inequality pi > p t for the pressures, contrary to the correct result pi < pt (for all relevant magnetic 
field strengths and densities, the Maxwell term is larger than the term linear on magnetization in 
Eq. (|18p . as was clarified earlier [3(3] and is confirmed by our calculations). Eqs. (fTTl) . (ITS]) were 
obtained in Ref. using the functional integration technique of a quantum field theory |39(|. Also, 
the similar equations (after disregarding the bag pressure and small higher order terms containing 



M) were obtained in Refs. [32 . 1331 ] for a degenerate system of strongly interacting neutrons in an 



external magnetic field by applying a Fermi liquid approach 4CJ, |41[. Moreover, the Maxwell term 
was also accounted for in the stress tensor of a relativistic magnetized system of free neutrons in 



Ref. 3l|, where a similar conclusion about the possibility of the collapse of a neutron star along 



the magnetic field (not the transverse collapse) was reached. 

Note that Eqs. (|16j) - ([18j) were obtained in Ref. |30j] for an abstract magnetized one-component 
degenerate Fermi gas. All calculations in that work were done at fixed chemical potential and 
zero bag pressure. We consider a scenario which is relevant in the astrophysical context of the 
problem, namely, we study the effects of the pressure anisotropy explicitly in magnetized three- 
flavor quark system (magnetized strange quark matter) under the additional constraints of the 
total baryon number conservation, charge neutrality and chemical equilibrium with respect to 
various weak processes occurring in the system. The quark species are characterized by their own 
chemical potentials and, together with the electron chemical potential, are determined from the 
constraints dH), ©, (| 13|) and (j!4p at the given total baryon number density. Under such formulation 
of the problem, it is possible to find how the chemical composition of strange quark matter, spin 
polarization of various quark flavors and EoS change with the magnetic field, putting a special 
emphasis on the effect of the pressure anisotropy on the above dependences. 



II. NUMERICAL RESULTS AND DISCUSSION 



As was mentioned in Introduction, strange quark matter can be in absolutely stable state 
(strange quark stars), or in metastable state, which can be stabilized by high enough external 
pressure (hybrid stars). The valley of the absolute stability in the model parameter space is 
determined from the requirement that, at zero external pressure and temperature, the energy per 
baryon for strange quark matter should be less than that for the most stable 56 Fe nucleus being 
about 930 MeV. The maximum allowed bag pressure from the absolute stability window decreases 
with the magnetic field strength [13] and reaches its maximum B ~ 90 MeV/fm 3 at H = 0. 
In turn, the lower bound on the bag pressure is established from the requirement that, at zero 
temperature and pressure, two-flavor quark matter (composed of u and d quarks) should be less 
stable with respect to the iron nucleus 56 Fe, and, hence, the energy per baryon for two-flavor quark 
matter should be larger than that for the nucleus 56 Fe Further, we will be interested in the 

astrophysical scenario, in which strange quark matter is formed in the core of a strongly magnetized 
neutron star, and, hence, is metastable at zero pressure. In numerical calculations, we adopt two 
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FIG. 1. (Color online) Various fermion species chemical potentials as functions of the magnetic field strength 
at zero temperature for the total baryon density (a) qb = 3po and (b) qb = 4go- The chemical potential 
of positrons is shown by the red curves between the vertical dotted lines. The vertical arrows indicate the 
points corresponding to the critical field H c ; see further details in the text. 



values of the bag constant, B = 100 MeV/fm 3 and B = 120 MeV/fm 3 , which are slightly larger 
than the upper bound on B from the absolute stability window. The core densities corresponding 
to these bag pressures are chosen equal to qb = 3go and qb = ^Qo, respectively, which are, in 
principle, sufficient to produce deconfinement [1] (go = 0.16 fm -3 being the nuclear saturation 



density). For the current quark masses, we use the same values as in Refs. [25l . |26| . |29| . |43J, i.e., 
m u = rrid = 5 MeV, and m s = 150 MeV. 

Fig. [TJ shows the chemical potentials of all fermion species as functions of the magnetic field 
strength. It is seen that the chemical potentials of fermions, first, stay practically constant under 
increasing the magnetic field, with d and s quark chemical potentials being somewhat larger (on 
the value of \x e - ~ 14 — 16 MeV) than the u quark chemical potential. The apparent Landau 
oscillations of the chemical potentials appear beginning from H ~ 3 ■ 10 18 — 4 • 10 18 G, depending 
on the total baryon number density. At H > 4 ■ 10 19 G, the quark chemical potentials decrease 
with the magnetic field. An interesting peculiarity occurs in a narrow interval near H ~ 2 ■ 10 19 G, 
marked by the vertical dotted lines. Namely, for magnetic field strengths from that interval the 
u quark chemical potential is larger than the d and s quark chemical potentials, \i u > = /i s . 
Hence, for such magnetic fields, according to Eq. (fl"3l) . the electron chemical potential would be 
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negative, fi e - < 0. If to recall the finite temperature expression for the electron number densitjQ, 
its zero temperature limit at \i e - < is, formally, negative, contrary to the constraint Q e - ^ 0. In 
fact, this means that in this interval on H electrons are missing and, hence, the weak j3~ processes 
(fTU|) . (fTTj) are impossible. However, for such magnetic fields, the following weak (3 + processes 
become allowable 

u^>-d + e + + v e , d + e + ^u + u e , (20) 
u — > s + e + + u e , s + e + — > u + D e . (21) 

Hence, for this specific range of the magnetic field strengths, the charge neutrality and chemical 
equilibrium conditions should read 

1q u ~ Qd ~ Qs + 3p e + = 0, (22) 
= fJ-d + Me+J fJ-d = f^s, (23) 
which should be solved jointly with the condition of the total baryon number conservation, Eq. ([H]). 
The quark and positron chemical potentials obtained as solutions of these equations are shown 
graphically in Fig. Q] as the corresponding curves between the vertical dotted lines. With increasing 
the core density, the width of the interval on H, where positrons appear, increases slightly as 
well (cf. the ranges 1.56 • 10 19 G-1.80 • 10 19 G at g B = 3g and 1.86 • 10 19 G-2.21 • 10 19 G at 
Qb — 4^0 )■ Thus, as a matter of principle, in strongly magnetized strange quark matter at zero 
temperature, subject to the total baryon number conservation, charge neutrality and chemical 
equilibrium conditions, positrons can appear in a certain narrow interval of the magnetic field 
strengths, replacing electrons. In this case, strange quark matter will have negative hadronic 
electric charge. Note that, according to Ref. [3j, the contact of stable strange quark matter, having 
negative hadronic electric charge, with the ordinary matter would have the disastrous consequences 
for the latter, because positively charged nuclei would be attracted to strange quark matter and 
absorbed. However, the contact of metastable strange quark matter, having negative hadronic 
electric charge, with hadronic matter in the interior of a neutron star is possible, because the outer 
hadronic layer provides the necessary external pressure to stabilize strange quark matter in the 
core and cannot be completely depleted. Nevertheless, we should calculate the critical field H c for 
the appearance of the longitudinal instability, which could prevent the occurrence of positrons in 
a certain range of magnetic field strengths with H > 10 19 G. The meaning of the vertical arrows 
in Fig. [I] will be discussed later in the text. 

Fig. [2] shows the abundances of various fermion species as functions of the magnetic field 
strength. The number densities of u and d quarks are quite close to each other for all mag- 
netic fields under consideration. The electron number density begins quite rapidly to increase at 
H « 2.2 • 10 16 G for g B = 3g and at H « 1.9 • 10 16 G for g B = Ag . As noted earlier, in the 
narrow interval near H ~ 2 • 10 19 G electrons are replaced by positrons, and beyond this inter- 
val electrons appear again with the number density increasing with H. The s quark content of 
strange quark matter stays practically constant till the field strength H ~ 4.1 • 10 18 G at g B = 3go 
and H ~ 3.8 • 10 18 G at g B = &g$, beyond which the s quark number density experiences visible 
Landau oscillations. Then, beginning from the field strength H ~ 3.2 • 10 19 G at g B = 3go and 
H ~ 3.9 • 10 19 G at g B = Ago, the s quark content rapidly decreases. Strange quark matter loses 
its strangeness and turns into two-flavor quark matter in the magnetic fields slightly larger than 

1 With account of the spin degeneracy factor (2 — 5u,o), it reads [25l |: 

„=n JO 



2tt 2 ^ ''Jo le^-^-'+l 



e /3( £ J+M c -) +1 
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10 20 G. Again, we should determine the critical field H c in order to check whether this significant 
drop of strangeness could happen in a strong magnetic field. 

Fig. [3] shows the spin polarization parameter IT for various fermion species, determined ac- 
cording to Eq. (jjj), as a function of the magnetic field strength. Spin polarization of u quarks is 
positive while for d, s quarks and electrons it is negative. The magnitude of the spin polarization 
parameter IT increases with H till it is saturated at the respective saturation field H^. At H = HI, 
the corresponding ith. fermion species becomes fully spin polarized. The respective values of the 
saturation field are: H e s » 2.2-10 16 G for electrons, H% « 1.7-10 19 G for u quarks, H s s w 2.5-10 19 G 
for s quarks and Hg ~ 3.2 • 10 19 G for d quarks at qb = 3£>0) an d H% ~ 1.9 • 10 16 G for electrons, 
H% w 2.0 • 10 19 G for u quarks, H s s « 3.1 • 10 19 G for s quarks and H% m 3.9 • 10 19 G for d quarks at 
Qb = 4^o- Note that quite a rapid increase of the electron number density with the magnetic field 
(cf. Fig. |2J) begins just at the saturation field Hf, and, hence, this increase occurs when electrons 
become completely spin polarized. Further oscillations in the electron number density are, in fact, 
caused by the Landau oscillations of the quark number densities, which influence the electron pop- 
ulation through the charge neutrality condition. The significant drop of the strange quark content 
begins just at the magnetic field strength at which the d quarks become completely spin polarized. 
Although the s-quark current mass is larger than that for d quark, m s > uid, s quarks become fully 
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FIG. 3. (Color online) Same as in Fig. [T]but for the spin polarization parameter of various fcrmion species. 
The spin polarization parameter of positrons is shown by the red segments between the vertical dotted lines. 



polarized at a smaller saturation field because their particle density is smaller than for d quarks, 
q s < Qd- The spin polarization parameter of various fermion species in the magnetic field range, 
where positrons appear, is shown by the respective curves between the vertical dotted lines. It is 
seen that positrons occur already fully polarized, and u quarks become totally polarized just in 
this range of the magnetic field strengths. Nevertheless, as mentioned before, only after determin- 
ing the critical field H c for the appearance of the longitudinal instability, it would be possible to 
determine the degree of spin polarization which could be reached for each of the fermion species. 

Now we present the results of calculations of the longitudinal pi and transverse pt pressures. 
Fig. H] presents these quantities as functions of the magnetic field strength at qb = 3qq,B = 
100MeV/fm 3 (solid lines) and qb = 4:Q ,B = 120MeV/fm 3 (dashed lines). It is seen that the 
transverse (ascending branches) and longitudinal (descending branches) pressures, first, stay prac- 
tically constant and indistinguishable from each other. This behavior of pt and p\ corresponds 
to the isotropic regime. Beyond some threshold magnetic field Hth, the transverse pressure 
Pt increases with H while the longitudinal pressure p\ decreases with it, clearly reflecting the 
anisotropic nature of the total pressure in strange quark matter in such strong magnetic fields 
(anisotropic regime). In the critical magnetic field H c , the longitudinal pressure p\ vanishes. This 
happens at H c « 7.4 • 10 17 G for q b = 3q ,B = 100MeV/fm 3 , and at H c « 1.4 • 10 18 G for 
qb = 4:Qo,B = 120MeV/fm 3 . Above the critical magnetic field, the longitudinal pressure is nega- 
tive leading to the longitudinal instability of strange quark matter. Therefore, the thermodynamic 
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H[G] 

FIG. 4. (Color online) Transverse (ascending branches in the branching curves) and longitudinal (descending 
branches in the branching curves) pressures in magnetized strange quark matter at zero temperature as 
functions of the magnetic field strength for qb = 3go,B = 100 MeV/fm 3 (solid lines) and qb = 4:Qo,B = 
120 MeV/fm 3 (dashed lines). 



properties of strange quark matter should be considered in the magnetic fields H < H c . In fact, 
the critical field sets the upper bound on the magnetic field strength which can be reached in the 
core of a strongly magnetized hybrid star. For comparison, we present also the values of the crit- 
ical field for dense neutron matter with the Skyrme BSk20 interaction at zero temperature being 
H c » 1.6 • 10 18 G at Q B = 3£> , and H c « 2.4 • 10 18 G at g B = 4g Q. 

Now, in order to see, which of the discussed already features of strange quark matter at zero 
temperature in a strong magnetic field are preserved before the appearance of the longitudinal 
instability, we show in Figs.[T][3]by the vertical arrows the respective values of the physical quantities 
corresponding to the critical field H c . Let us begin with Fig.[T]for the chemical potentials of various 
fermion species. It is seen that the chemical potentials of quarks and electrons stay practically 
unchanged before the appearance of the longitudinal instability. The significant changes in the 
chemical potentials occur only in the fields H > H c . In particular, the longitudinal instability 
precludes the appearance of positrons for which the fields H > 10 19 G are necessary. 

Let us turn to Fig. [2] for the abundances of various fermion species. Till the critical field H c , the 
content of quark species stays practically constant while the electron fraction remains quite small, 
Qe~ 1 00 10~ 2 . Also, there is no room for the significant drop of the strange quark content in strong 
magnetic fields H ~ 10 20 G which is averted by the appearance of the longitudinal instability in 
the critical field H c . In fact, despite the presence of strong magnetic fields H ~ 10 18 G, strange 
quark matter has the same fraction of s quarks as in the field-free case. 

Let us now consider Fig. [3] for spin polarizations of various fermion species. It is seen that the 
full polarization in a strong magnetic field can be achieved only for electrons. For various quark 
species, the spin polarization remains quite moderate up to the critical magnetic field H c . E.g., at 
qb = ^Qo,H = H c we have IL U sa 0.06, ILi pa —0.03, II S « —0.04; at qb = 3qq,H = H c , the quark 
spin polarizations are similar to these values with the maximum magnitude of the spin polarization 
parameter for u quarks, ILj ~ 0.04. Therefore, the occurrence of a field-induced fully polarized 
state in strange quark matter is prevented by the appearance of the longitudinal instability in the 
critical magnetic field. The degree of spin polarization of various constituents is an important issue 
for determining the neutrino mean free paths in magnetized strange quark matter 44j |. and, hence, 
it is relevant for the adequate description of the neutrino transport and thermal evolution of a 
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FIG. 5. (Color online) Same as in Fig. |4]but for the normalized difference 6 = Pt ~f l between the trans- 
verse and longitudinal pressures. The vertical arrows show the maximum normalized splitting 5 C at the 
corresponding critical field H c . 



magnetized hybrid star. 

Thus, in the anisotropic regime, the pressure anisotropy in a strong magnetic field plays an 
important role and should be accounted for in the description of the thermodynamic properties of 
strange quark matter. Let us now make the estimate of the threshold magnetic field Hth above 
which the pressure anisotropy cannot be disregarded. Fig. [5] shows the normalized difference 
between the transverse and longitudinal pressures 

x Pt-Pl 

o = , 

Po 

where po is the isotropic pressure (which corresponds to the weak field limit with pi = pt = pp). as a 
function of the magnetic field strength for the cases under consideration. Following Refs. 



for finding the threshold field H t h one can use the approximate criterion 6 ~ 1. Then anisotropic 
regime enters at H th « 5.5 • 10 17 G for q b = 3q ,B = 100MeV/fm 3 , and at H th m 9.9 • 10 17 G 



for qb = 4^0 , B = 120MeV/fm 3 . For comparison 33j, the threshold field for neutron matter 
with the BSk20 Skyrme interaction at zero temperature is H^ ~ 1.2 • 10 18 G for qb = 3qo, and 
H t h ~ 1.8 • 10 18 G for qb = 4^o- The anisotropy parameter 5 reaches its maximum S c ~ 2 in 
the critical field H c , corresponding to the onset of the longitudinal instability in strange quark 
matter. In the anisotropic regime, a hybrid star is deformed and takes the oblate form. Thus, 
as follows from the previous discussions, the effects of the pressure anisotropy are important at 
Hth < H < H c , and significantly influence the structural and polarization properties of the quark 
core in a strongly magnetized neutron star. 

Fig. [6] shows the energy density E of the system without the magnetic field energy density con- 
tribution Ef = (top panel) and with account of it (bottom panel) as a function of the magnetic 
field strength at zero temperature. It is seen that, due to the Landau diamagnetism, the energy 
density of solely magnetized strange quark matter decreases with the magnetic field. However, the 
overall effect of the magnetic field, with account of the pure magnetic field contribution Ef, is to 
increase the energy density of the system. Nevertheless, this effect of the magnetic field is, in fact, 
insignificant because the magnetic field is bound from above by the critical magnetic field H c . The 
values of the energy density E, corresponding to the critical field H c , are shown in Fig. [6] by the 



12 



1200 



> 



LU 

I 

LU 



1000 
800 
600 
400 
200 



P B = 3p , B=100 MeV/fm 

p =4p n , B=120 MeV/fm 3 





10 1 



1400 
1200 
1000 



> 800 



LU 600 
400 



200 



(a) 



10' 



10 18 



H[G] 



10' 



10' 



10 18 



H[G] 



10 



'_ p B = 3p , B=100 MeV/fm 3 


1 f 
1 I 
t 1 
1 f 


■ p B = 4p , B=120 MeV/fm 3 

J. 


I 1 
I I 
i f 
I I 
I I 
/ I 
t I 


±-^- 




■ (b) 

' ' 





10' 



FIG. 6. (Color online) Same as in Fig. [4] but for the energy density E of the system (a) without the magnetic 
field energy density contribution Ef = ^ and (b) with account of it. The vertical arrows indicate the points 
corresponding to the critical field H c . 



vertical arrows. 

Because of the pressure anisotropy, the equation of state of strange quark matter in a strong 
magnetic field is also anisotropic. Fig. [7]shows the dependence of the energy density E of the system 
on the transverse pressure pt (top panel) and on the longitudinal pressure pi (bottom panel) after 
excluding the dependence on H in these quantities. In particular, the anisotropic character of 
the pressure is reflected in the fact that the energy density is the increasing function of pt while 
it decreases with pi. This is because the dominant Maxwell term enters the transverse pressure 
Pt and the energy density E with positive sign while it enters the longitudinal pressure pi with 
negative sign. The vertical arrows in the top panel indicate the points in these lines corresponding 
to the critical field H c . In the bottom panel, the physical region corresponds to pi > 0. 

Note that because the EoS of strange quark matter becomes essentially anisotropic in an ul- 
trastrong magnetic field, the usual scheme for finding the mass-radius relationship based on the 
Tolman-Oppenheimer-Volkoff (TOV) equations |45(] for a spherically symmetric and static com- 
pact star should be revised. Instead, the corresponding relationship should be found by the self- 
consistent treatment of the anisotropic EoS and axisymmetric TOV equations substituting the 
conventional TOV equations in the case of an axisymmetric compact star. 
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FIG. 7. (Color online) The energy density E of the system at zero temperature as a function of: (a) the 
transverse pressure p t and (b) the longitudinal pressure for the cases qb = 3qo,B = 100 MeV/fm 3 and 
qb = 4go, B = 120 MeV/fm 3 . The meaning of the vertical arrows in the top panel is the same as in Fig. [6] 
In the bottom panel, the physical region corresponds to pi > 0. 



In summary, we have considered the impact of strong magnetic fields up to 10 G on the ther- 
modynamic properties of strange quark matter at zero temperature under additional constraints 
of total baryon number conservation, charge neutrality and chemical equilibrium with respect to 
various weak processes occurring in the system. The study has been done within the framework 
of the MIT bag model with the finite current quark masses m u = ^ m s . In the numeri- 
cal calculations, we have adopted two sets of the total baryon number density and bag pressure, 
Qb = 3qq,B = 100 MeV/fm 3 and qb = ^Qo,B = 120 MeV/fm 3 . It has been found that in strong 
magnetic fields up to 10 20 G some interesting features in the chemical composition and spin struc- 
ture of strange quark matter could occur: 

(1) The content of strange quarks rapidly decreases in the fields somewhat larger than 10 19 G 
and becomes negligible in the fields slightly exceeding 10 20 G; 

(2) For the magnetic field strengths in the quite narrow interval near H ~ 2 ■ 10 19 G the 
constraints of total baryon number conservation, charge neutrality and chemical equilibrium can 
be satisfied only if positrons appear in various weak processes in that range of the field strengths 
(instead of electrons); 
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(3) Electrons occupy only the lowest Landau level and, hence, become completely spin polarized 
in the magnetic fields somewhat larger than 10 16 G; u, s and d quarks become fully polarized in 
the fields somewhat larger than 10 19 G (the recitation of the quark species is in the order in which 
they appear fully polarized under increasing H). 

Nevertheless, under such strong magnetic fields, the total pressure containing also the magnetic 
field contribution, becomes anisotropic, and the effects of the pressure anisotropy change most 
of the above conclusions. Namely, the longitudinal (along the magnetic field) pressure decreases 
with the magnetic field (contrary to the transverse pressure increasing with H) and vanishes in 
the critical field H c resulting in the longitudinal instability of strange quark matter. The value 
of the critical field H c depends on the total baryon number density of strange quark matter and 
the bag pressure B, and it turns out to be somewhat less or larger than 10 18 G for the two sets 
of the parameters, considered in the given study. Therefore, the appearance of the longitudinal 
instability in strong magnetic fields beyond the critical one precludes the features (1), (2) in the 
chemical composition of strongly magnetized strange quark matter. Concerning the conclusion 
(3), only electrons can reach the state of full polarization, that is not true for quarks of all flavors, 
whose polarization remains mild even for magnetic fields near H c . 

The pressure anisotropy becomes relevant beyond some threshold field H t h- For the sets of the 
total baryon number density and bag constant considered in the given study, it turns out that 
10 17 < Hfh ^$ 10 18 G. This estimate is somewhat less than that found for strongly magnetized 
dense neutron matter, H t h ~ 10 18 G [H]. In strong magnetic fields H > Hth, the EoS of 
strange quark matter becomes essentially anisotropic. The longitudinal and transverse pressures 
as well as the anisotropic EoS of magnetized strange quark matter have been determined at the 
total baryon number densities and magnetic field strengths relevant to the interiors of magnetars. 

In this work, we have studied the impact of a strong magnetic field on the thermodynamic 
properties of strange quark matter at zero temperature. It would be also of interest to extend this 
research to finite temperatures 33, |4fj], which can lead to a number of interesting effects, such as, 
e.g, an unusual behavior of the entropy of a spin polarized state 47j, 1481 ]. 

In conclusion, it is worthy to note that, because strong mag netic fields of about H ~ 10 18 G 
(RHIC), or even by order of magnitude larger (LHC), are generated in non-central high-energy 
heavy- ion collisions [49l.l50l|. the effects of the pressure anisotropy should be relevant there as well. In 
particular, the pressure anisotropy in a strong magnetic field can contribute to the enhancement of 
the elliptic flow of hot nuclear matter created in a heavy-ion collision 5l| . Because the conditions in 
high-energy heavy-ion collisions are different from those in the cores of strongly magnetized neutron 
stars (the absence of chemical equilibrium with respect to the weak processes, the low baryon 
number densities and high temperatures), the possibility for the occurrence of the strangeness 
suppression in a strong magnetic field for the former case needs a separate study. 

J.Y. was supported by grant 2010-0011378 from Basic Science Research Program through NRF 
of Korea funded by MEST and by grant R32-10130 from WCU project of MEST and NRF. 
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